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Abstract 

We present a semiclassical analysis of the instability of an electron shuttle composed of three quantum dots: two are fixed 
and coupled via leads to electron resevoirs at u.^ r with u.^ S> (j.r, while the central dot is mounted on a classical harmonic 
oscillator. The semiclassical analysis, which is valid if the central dot oscillation amplitude is larger than the quantum 
mechanical zero point motion, can be used to gain additional insight about the relationship of resonances and instabilities 
of the device. 
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In nano-electromechanical systems (NEMS) the elec- 
trical and mechanical properties are deeply intercon- 
nected. An archetypal NEMS device consists of a mov- 
able object connected to leads [1] . The charge distribu- 
tion gives rise to an electrical force, which influences the 
mechanical dynamics, while the position of the mov- 
able object determines the tunneling rates from the 
leads and thus influences the electrical dynamics of the 
system. A current through the device can sustain me- 
chanical oscillations even in the presence of damping. 
An interesting regime of transport arises when only 
one electron per cycle is transferred from the left to 
the right lead. Due to the position dependent tunneling 
amplitude the movable part gets charged when near to 
the left lead, then the electrostatic force pushes it to- 
wards the right lead where the now enhanced tunneling 
rate helps the release of the electron. 

We describe the electronic part with the density ma- 
trix formalism and couple the master equation to a 



classical equation of motion for the central dot posi- 
tion [1]. We perform a linear instability analysis of the 
system to pinpoint where the equilibrium solution for 
the system becomes unstable leading to Hopf bifurca- 
tions; the relative phase between charge, velocity and 
position of the unstable rotating solution shows the 
features of shuttling in agreement with the quantum- 
phase-space description [2] . 

The semiclassical approach is justified by the 
quantum-classical correspondence since the dot oscil- 
lations are bigger than the minimum quantum ampli- 
tude. However the mean field approach for the electric 
part neglects the effect of shot noise and the damping 
factor threshold is usually very much reduced com- 
pared to the quantum treatment. The semiclassical 
analysis is also much easier to handle numerically. 
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1. The model 

Consider a simple NEMS consisting a small metallic 
grain (or a QD) connected to two leads. The grain is 
moving in a parabolic potential. Two additional dots, 
one on each side of the oscillating component, fix the 
energy of the incoming and outgoing electrons [3] . The 
Hamiltonian of the system is H — Hm ec h + H e i , where 



H — H„ 



where 



Hmech 
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We assume strong Coulomb blockade regime, thus the 
vectors = L, C, R, together with the empty state 
|0), span the entire single particle Hilbert space of the 
device. The matrix elements are explicitly: 
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AV, the device bias, is the difference between the en- 
ergy of the left and the right dot. xo is half the dis- 
tance between the two outer dots and represents the 
maximum amplitude of the inner dot oscillation. The 
three dots are electrically connected only via a tunnel- 
ing mechanism. The tunneling length is given by 1/a 
and the tunneling strength depends on the position x 
of the inner grain, 

t L (x) = V e- a(xo+x) 

i i ( 4 ) 
t R (x) = Voe- a(x °- x) 

The electric part of the NEMS is described in the den- 
sity matrix formalism. The Hilbert space has dimen- 
sion four, though the corresponding four-by-four den- 
sity matrix can be reduced to an effective three-by- 
three matrix since the elements of the form po; , i = 
L,C,R are decoupled from the others and poo can be 
eliminated using J2iP" = 1- This effective density ma- 
trix obeys the generalized master equation: 



p = -i[H,p] + Hp 



(5) 



The first term represents the coherent evolution of the 
electrons in the three dots when isolated from the leads. 



The coupling to the leads is introduced in the wide 
band approximation following Gurvitz [4] and gives the 
second term of the equation (Hp): 
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where V is the injection rate to the leads. The equa- 
tion of motion for the central dot position has a con- 
servative part deduced from the (quantum-) averaged 
Hamiltonian (H) — H mech + Tr(pH e i(x)) and a phe- 
nomenological damping term: 

^_d{H) = p_ 
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2. Instability analysis 

In the shuttling regime the electrical and the me- 
chanical component of the system evolve with the same 
time scale. Thus it is important to treat all degrees of 
freedom on equal footing: we introduce a set of gener- 
alized coordinates q, and write the Eqs. (5,7) as 

Qi = fi(q) (8) 

and linearize around the equilibrium, (0 = F(q )): 

qi = Mij(qj - q°) + ... (9) 



where Mij = djFi\ q=q 



In general, even if the veloc- 
ity of the central dot is zero, a current can flow through 
the system and is given by I = TpuR. Let us study 
the equilibrium position and electronic configuration 
of the three dot system as a function of the device bias. 
At AV = the central dot is at rest in the middle (we 
will take this as zero position). This unique solution 
appears natural, but in fact is not obvious if one rec- 
ognizes the presence of an exchange force that acts on 
the dot (i.e. the last line of Eq. 7) 2 . As the difference 
between the energy of the left and right dot is increased 
the equilibrium position of the central dot shifts to the 
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right, but then returns to zero again for high device 
biases (Fig. la). This reveals a competition between 
the electrostatic force and the progressively dominat- 
ing decoupling of the three dots due to the mismatch 
of their energy levels. The equilibrium occupation of 
the three dots also reflects this decoupling (Fig. lb). 
In the limit of high device biases the left dot has the 
highest probability of being occupied while the others 
are almost empty. At small biases the equilibrium oc- 
cupation is strongly dependent on the injection rate V. 
In this regime the ratio between the injecting rate and 
the bare tunneling probability r — T/(Voe~ ax °) accu- 
rately describes the scenario. If r -C 1 the three dots 
are almost isolated and their coherent oscillations are 
damped by the contacts, but so slowly that the three 
dots tend to a uniform stationary occupation. In the 
opposite limit (r S> 1) the left dot is continuously re- 
filled and the right emptied at a rate which overcomes 
the dot's dynamical response. The left and central dot 
share the occupation probability while the right dot 
is empty and the empty state occupation probability 
vanishes. 




Fig. 1. Equilibrium position (a) and equilibrium electronic 
configuration (b) as a function of the device bias for three 
different injection rates T = 0.05,0.5,5. 

The quantum description of the shuttle device shows 
that the central dot can oscillate with an amplitude 
which is much bigger than the minimum uncertainty 
length [2]. We expect the equilibrium solution of the 
semiclassical description to be unstable in the same 
parameter range (except for the damping: the lack 



of noise sources in the present semiclassical approach 
reduces the instability threshold). The vector q = 

[ pLL PCC PRR PLC PCL pLR pRL PCR PRC X P } has 

an electrical part (the first 9 components) and a me- 
chanical part (the last 2 components). The matrix M 
is written in a block form, consistent with this natural 
separation: 



M 



A C 
V B 



(10) 



The eigenvalues of M represent the characteristic fre- 
quencies and damping for small oscillation of the po- 
sition and the density matrix around the equilibrium 
configuration. If at least one of the eigenvalues of M 
has a positive real part then the equilibrium solution 
is unstable. In the general case all the four blocks of 
the matrix M depend on both the electrical and me- 
chanical degree of freedom. But if we set xq = 1/a and 
take the limit a —> 0, the off-diagonal blocks of the M 
(C and T>) vanish and the eigenvalues are those of A 
(electrical) and B (mechanical) . The imaginary part of 
the spectrum consists of the Bohr frequencies (i.e. all 
possible differences between the eigenvalues) of the iso- 
lated electron system and the mechanical frequency of 
the oscillator (Fig.2a). The real part of the spectrum, 
though, is negative due to the presence of the leads and 
the mechanical damping (Fig. 2b). As a is increased, 
the spectrum of A4 gets a positive real part: two in- 
stability regions appear (Fig. 2b). They have the form 
of Hopf bifurcations: for these parameters the solution 
of Eq.8 is a limit cycle. Correspondingly, in the unsta- 
ble regime, the imaginary part of the spectrum shows 
crossings between the "electrical" and "mechanical" 
frequencies (arrows in Fig. 2a). These are resonances 
between the oscillator and the electrical modes of the 
three dots. They appear AV ~ ui and 2lo. The slight 
shift from the exact values is due to the bare interac- 
tion Vb which modifies the electronic spectrum of the 
Hamiltonian from the decoupled 0, ±/\V/2 (which re- 
mains as asymptotic behavior). 

Only an analysis of the eigenvectors can demon- 
strate the actual oscillation of both the electric and 
mechanical components in the limit cycle solution. For 
this reason we have studied the electrical and mechan- 
ical weight of the eigenvectors in the unstable regime 
(Fig.3a). The larger the real part of the unstable eigen- 
value, the higher is the mixing of the two components. 
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In the quantum description [2] , the shuttling insta- 
bility is characterized by a strong correlation between 
charge, position and velocity of the oscillating dot. 
Therefore we study next the relative phase between 
charge, velocity and position of the unstable solution, 
even if a priori is it not clear whether this property will 
be maintained in the limit cycle stationary solution. 
The phase analysis shows, for a given charge config- 
uration (that for definiteness we take maximal) a de- 
crease in the position and velocity phase: for both res- 
onances the charged dot passes from maximal position 
and zero velocity to minimal position and zero veloc- 
ity. The range of phase rotation is decreasing together 
with the instability in the case with A V ~ uo (Fig. 3b) . 
This rotation of the charged dot configuration in phase 
space is also in qualitative agreement with the quan- 
tum analysis. 
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Fig. 3. Eigenvector analysis as a function of device bias, a) 
The relative amplitude of the electrical (full) and mechan- 
ical (dot dashed) components of the unstable eigenvector. 
The relative amplitude of the two components tends to be 
equal in the instability region and the effect is more promi- 
nent for higher a. b) The phase of position (full) and ve- 
locity (dot dashed) with respect to the phase of the charge 
in the central dot for the unstable eigenvector. 

framework to identify the correspondence between in- 
stability and resonance in the NEMS. Also the rota- 
tion of the phase space configuration of the charged 
dot in the shuttling regime was detected. More work is 
needed to understand whether other characteristics of 
the system such the higher order instabilities (AV « 
3u},4ui, ...) are pure quantum phenomena. 
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Fig. 2. Imaginary (a) and real (b) part of the spectrum of A4 
for different values of tunneling length, a = 0,0.1,0.2,0.5. 
The instability grows with a. Only the four eigenvalues 
with positive imaginary part are plotted. The system is 
symmetric to complex conjugation. The three eigenvalues 
with zero imaginary part are omitted for simplicity: they 
represent always damped modes. 
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3. Conclusions 



In this work we have analyzed the instabilities of a 
triple-dot shuttle using a semiclassical approach, moti- 
vated by the "semiclassical behavior" of the quantum 
description. We were able in this simple and intuitive 
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